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Abstract 


The present report describes the development of an interaction data base and a numerical 
solution to the transport of baryons through an arbitrary shield material based on a straight 
ahead approximation of the Boltzmann equation. The code is most accurate for continuous 
energy boundary values but gives reasonable results for discrete spectra at the boundary with 
even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H 2 O). 

1. Introduction 

The purpose of the present report is to describe computer programs developed for the 
calculation of the transport of high energy nucleons (baryons) and their interaction with 
materials. The methods, based on the direct solution of the Boltzmann equation, have been 
developed over the last several years (refs. 1-4). The solutions employ the straight ahead 
approximations for which energy changing processes are accurately handled, but all propagation 
occurs along a fixed direction. The present goal is to document a relatively complete description 
of the basic physical processes even though the input data base used requires some improvement. 
Future work will concentrate on improving the data base, code efficiency, and the computational 
procedures. 


2 . Theoretical Considerations 


2.1. Boltzmann Equation 

The equations in the straight ahead approximation to be solved (ref. 5) are 






__ roc 
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where <j>j(x,E) is the type j particle differential flux density at x with energy E\ S(E) is 
proton stopping power; a p (E ), a n (E) are proton, neutron total cross section, respectively; 
and fij{E, E f ) are the differential cross sections for elastic and nonelastic processes. New field 
quantities are conveniently defined as 
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which can be rewritten as integral equations with boundary at x = 0. The results follow: 


ipp{x, r) = e fo a p ( r+z ) dz ij} p (0, r + x) 
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These equations relate the known flux at the boundary at x = 0 to an 
choose the point to be x = h , where h is a small distance, and write 
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(2.1.11) 

To evaluate equations (2.1.10) and (2.1.11), the flux ipj(x, r) is required 
This we do as follows: 

If h is sufficiently small such that 

across the interval. 

Cj(r')h < 1 

(2.1.12) 

then according to perturbation theory (ref. 1) 
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which may be used to approximate the above integrals of equations (2.1.10) and (2.1.11). 

For many cases of practical interest (e.g., accelerator studies), monoenergetic particle beams 
are used, and separation of the singular terms from the solution becomes convenient. The initial 
beam of type J particles of energy £o(ro = R{Eq)) is taken as 

ipj(0,r) = 6jj6(ro - r ) 

(2.1.15) 

and the solution is written as 


rpj(x, r) = tl>j 0 (x, r) + tpj(x, r) 

(2.1.16) 

The corresponding singular terms are 
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The regular terms of equations (2.1.10) and (2.1.11) may be written as 
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The singular contributions under the integrals of equations (2.1.19) and (2.1.20) can be 
evaluated with equations (2.1.17) and (2.1.18) and the approximations (2. 1 . 12)— (2. 1 . 14) applied 
to find 
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where Tq = tq - x - § and F is related to the cumulative spectrum F as 
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r, r') = J /ij (r + z, r') dz 

= *ij(r + r') - F^(r,r') (2.1.23) 
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and e(r) is the energy associated with residual range r, and E f = e(r'). Equations (2.1.20) 
and (2.1.21) are evaluated by establishing an x-grid at which t/>y(z m ,r) is evaluated, and h 
is the distance between each successive evaluation. The integral over r f is accomplished by 
establishing an r - grid (and the corresponding E-grid) and using 


roo 

/ g{rn,r')ipj{x m ,r') dr' 

Jr n 


™ _ r r t + 1 
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(2.1.25) 


where fp = [rp + r^ + j)/2, and the series terminates at the highest t value in the r-grid. 
There is a spatially dependent discontinuity in the proton flux spectrum that requires right- 
and left-hand interpolation and integration. These discontinuities have been treated in the 
computational procedures. 


2.2. Neutron Source 

The neutron transport equation in three dimensions is 

poo 

[n- v + ( 7 n(E)]^„(f,n,E) = V / f nj {E,E\?i,Q!) do! dE' (2.2.1) 

j Je 

While it is true that the straight ahead approximation is adequate for most proton calculations, 
the neutron fields are more strongly affected by non-forward scattering components, particularly 
the low energy neutrons. The reason is that the lower energy neutrons have a greater range than 
lower energy protons because of the electric charge difference. Thus, a first-order correction 
to the straight ahead approximation may be applied by substituting the proton coupling in 
equation (2.2.1) by the straight ahead solution for 0 p (x, H,E) so that 

w f°° - - 

M 2 , n, e) = / f n p(E, e n, n x ) <t> p (x, e') dE' (2.2.2) 

Je 

The corresponding neutron transport equation is 

poo 

[n-V + e n {E)]<t>n(x,d,E) = / fnn(E,E',n,n') 4>n{x,Cl',E') dtl' dE' + tn(x,Q,E) (2.2.3) 

Je 

The neutron source integral is evaluated in a fashion similar to equation (2.1.25). 


2.3. Target Fragment Secondary Flux 

The target fragmentations produced in nuclear collision with the nucleon field must now be 
treated. The spectral parameters of the composite fragments are relatively independent of the 
projectile charge, energy, or direction. This leads to some simplifying assumptions so that 


d>: 




Et = RJ l 


Rj{E) + d(Q) 


(2.3.1) 

(2.3.2) 


where d(u) is the distance from the boundary (ref. 2). The source of ions of type j is evaluated 


as 
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where (Jji{E) are the fragmentation cross sections. The term fj{E) represents the spectral 
contributions averaged over all the target atomic constituents. In the present code, the distance 
to the boundary is assumed large. One could treat not only boundary effects but interface effects 
as well. 


3. Transport Coefficients 
3.1 • Stopping Power 

In passing through a material, an ion loses the larger fraction of its energy to electronic 
excitation of the material. Although a satisfactory theory of high-energy ion-electron inter- 
action is available in the form of Bethe’s theory utilizing the Born approximation, an equally 
satisfactory theory for low energies is not available. Bethe’s high-energy approximation to the 
energy loss per unit path (that is, stopping power) is given as 


S e = 


47rArZ 1 2 £ 2 e 4 

mv 2 



2mv 2 

{l~0 2 )h 



(3-1.1) 


where Z\ is the projectile charge, N is the number of target molecules per rmit volume, Z 2 is 
the number of electrons per target molecule, m is the electron mass, v is the projectile velocity, 
/? = v/c, c is the velocity of light, C is the velocity-dependent shell correction term (ref. 6), 
and I 2 is the mean excitation energy given by 


Z 2 Hh)=J2fnHEn) 

n 


(3.1.2) 


where f n are the electric dipole oscillator strengths of the target and E n are the corresponding 
excitation energies. Note that the sum in equation (3.1.2) includes discrete and continuum 
levels. Molecular stopping power is reasonably approximated by the sum of the corresponding 
empirically derived “atomic” stopping powers for which equations (3.1.1) and (3.1.2) imply 
that 

Z ln(7) = Ys n i z jWi) (3-1-3) 

3 

where Z and I pertain to the molecule, Zj and Ij are the corresponding atomic values, and rij 
are the stoichiometric coefficients. This additive rule (eq. (3.1.3)) is usually called Bragg’s rule 
(ref. 7). 

Sources of deviations from Bragg’s additive rule for molecules and for the condensed phase 
are discussed by Platzman (ref. 8). Aside from shifts in excitation energies and adjustments 
in line strengths as a result ,of molecular bonding, new terms in the stopping power appear 
due to coupling between vibrational and rotational modes. Additionally, in the condensed 
phase, some discrete transitions are moved into the continuum, and collective modes among 
valence electrons in adjacent atoms produce new terms in the absorption spectrum to be dealt 
with. Platzman proposed that the experimentally observed additivity rule may not show 
that molecular stopping power is the sum of atomic processes, but rather demonstrates that 
molecular bond shifts for covalent bonded molecules are relatively independent of the molecular 
combination. On the basis of such arguments, Platzman suggested that ionic bonded substances 
should be studied as a rigid test of the additive rule because of the radical difference in bonding 
type. He further estimated that ionic bond shifts could change the stopping power by as much 
as 50 percent. Recent results on shifts of mean excitation energies due to molecular bonds are 
discussed in references 9 to 11. Effects of the physical state have likewise been studied (ref. 12). 

The electronic stopping power for protons is adequately described by equation (3.1.1) 
for energies above 500 keV for which the shell or “tight binding” correction C makes an 
important contribution below 10 MeV (ref. 13). For proton energies below 500 keV, charge 
exchange (electron transfer) reactions alter the proton charge over much of its path, so that 
equation (3.1.1) is to be understood in terms of an average over the proton charge states. 


5 



Normally an average over the charge states is introduced into equation (3.1.1) so that the 
effective charge is the root-mean-square ion charge and not the average ion charge. At any 
ion energy, charge equilibrium is established very quickly in all materials. Utilizing the 
effective charge in equation (3.1.1) appears to make only modest improvement below 500 keV, 
an indication presumably of the failure of this theory based on an empirical basis (refs. 13 
and 14). The resultant stopping power for protons in water is shown in comparison with the 
evaluated data of Bichsel (ref. 15) in figure 1. 

The electronic stopping power for alpha particles requires terms in equation (3.1.1) of higher 
order in the projectile charge Z\ because of corrections to the Born approximation. The alpha 
stopping power cannot be related to the proton stopping power through their effective charges. 
Parametric fits to experimental data are given by Ziegler in reference 16 for all elements in 
both the gaseous and the condensed phases. 

The electronic stopping powers for heavier ions are related to the alpha stopping power 
through their corresponding effective charges. The effective charge suggested by Barkas (ref. 17) 
is used: 

Z* = z[l-exp(-125/?/£ 2/3 )] (3.1.4) 

where Z in equation (3.1.4) is the atomic number of the ion. 

At sufficiently low energies, the energy lost by an ion in a nuclear collision becomes 
important. The nuclear stopping theory used in this paper is a modification of the theory 
of Lindhard, Scharff, and Schiott (ref. 18). The reduced energy is given as 


e = 


32.53AiA 2 £ 


Z 1 Z 2 (A l +A 2 ) ( z ^ 3 + z | /3 ) 1/2 


(3.1.5) 


where E is in units of keV/nucleon, and A\ and A 2 are the atomic masses of the projectile and 
the target. The nuclear stopping power in reduced units . (ref. 16) is 


( 1.59c 1 / 2 (e < 0.01) ) 

* - < (0.01 < e < 10) 

(e > 10) J 

and the conversion factor to units of eV/10 15 atoms/cm 2 is 

8.426ZiA 2 Ai 


ln(0.47£) 

2t 


/ = 


(A l+ A 2 ) ( zj / 3 + z 2 2/3 ) 1/2 


(3.1.6) 


(3.1.7) 


The total stopping power Sj is obtained by summing the electronic and nuclear contribu- 
tions. Other processes of energy transfer such as Bremsstrahlung and pair production are 
unimportant. 

For energies above a few MeV/nucleon, Bethe’s equation is adequate provided that appro- 
priate corrections to Bragg’s rule (refs. 9 to 11), shell corrections (refs. 6, 13, and 14), and an 
effective charge are included. Electronic stopping power for protons is calculated from the para- 
metric formulas of Andersen and Ziegler (ref. 13). The calculated stopping power for protons 
in water is shown in figure 1 in comparison with data given by Bichsel (ref. 15). 

Because alpha stopping power is not derivable from the proton stopping power formula 
by using the effective charge at low energy, the parametric fits to empirical alpha stopping 
powers given by Ziegler (ref. 16) are used. Applying his results for condensed phase water 
poorly represented the data of references 19 and 20. Considering that physical state and 
molecular binding effects are most important for hydrogen (ref. 9), the water stopping power 
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was approximated by using the condensed phase parameters for hydrogen and the gas phase 
parameters for oxygen (which are known experimentally). These results are compared with 
experimental data for condensed phase water (refs. 19 and 20) in figure 1. It appears that 
Ziegler overestimated the condensed phase effects for oxygen since the gas phase oxygen data 
give satisfactory results as seen in figure 1. 

Electronic stopping powers for ions with a charge greater than 2 are related to the alpha 
stopping power through the effective charge given by equation (3.1.4). For water, the condensed 
phase formula of Ziegler for alpha particles gives probably the best stopping powers for heavier 
ions. Calculated results for 16 0 and 56 Fe ions in water are shown in figure 1 in comparison with 
the Northcliffe and Schilling data (ref. 21). Good agreement with Northcliffe and Schilling for 
56 Fe ions is especially important since their data seem to agree with the range experiments of 
J. H. Chan in Lexan resin (ref. 22). The stopping powers in Lexan resin and tissue equivalent 
material can be calculated in a way similar to the procedure given above for calculating stopping 
powers in water. 

3.2. Nuclear Absorption Cross Sections 

The nuclear absorption cross sections are calculated using a parameterization obtained by 
fitting an analytical expression to experimental data. For nucleon-nucleus collisions, it is given 
in millibarns by (ref. 23) 

a A = /(£)45A°' 7 {1 + 0.016 sin[5.3 - 2.63 ln(A)]} (3.2.1) 

where A is the mass number of the target nucleus, and the energy-dependent factor is 

f{E) = 1 - 0.62e~ E/20 ° sin (10.9-ET 028 ) (3.2.2) 


with the incident kinetic energy E in units of MeV/nucleon. 


3.3. Total Nuclear Cross Sections 

The proton-proton ( pp ) total cross section is found to be approximately 


<rpp{E) = (1 + 5/E ) 
for E > 25 MeV and taken as 


40+109 cos(0.199 


x/®)] 


,-.451(E-25)- 258 


_ ( F \ _ „ 6 . 51e-( E / 134 )' 7 

Vpp{E) - e 

at lower energies. The neutron-proton ( np ) cross section is taken as 

a np {E) = 38 + 12500e -1187 ( £-1 )' 35 


(3.3.1) 


(3.3.2) 


(3.3.3) 


which is used for E > 0.1 MeV and taken as 

a np (E) = 26000e _(£; /- 282) ' 3 (3.3.4) 


at lower energies. 

The low energy neutron-nucleus cross sections exhibit a complicated resonance structure 
over a nearly constant background. Above 20 MeV, the cross sections decrease with energy to 
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a minimum at 300 MeV and rise slightly to 1 GeV and remain nearly constant thereafter. In 
the present code, the total neutron cross sections are approximated at four energies as 


<r t (E,A)={ 


( 0.3A- 656 + 0.5.4e- ao66A 
0.38A 525 + 0.09Ae-° 5A 
0.45.4- 475 + 0.025j4e --06j4 
l 0.052.4- 726 + 0.02yle“ 025A 


(E = 0.1 MeV) ) 
(E = 1.5 MeV) 
(E = 20 MeV) 
(E = 100 MeV) J 

These are extended to lower and higher energies and interpolated using 

a(E,A) = <r(E i ,A)e~ a ^ E ~ E ^ 


(3.3.5) 


(3.3.6) 


where a (E{, A) is the appropriate cross section at one of the above energies and a is determined 
to ensure appropriate interpolation or extrapolation. This quantity is used to calculate the 
scattering cross section as 

o s {E) = amt {E) - <r abs (£) (3.3.7) 

used for crtot(£’) > 2o abs (£ : ), which occurs for values of E less than or about 200 MeV and 


cr s(-E') — CT abs(-^') 


(3.3.8) 


otherwise. 

The total neutron-nucleus cross section is shown in comparison with experimental data 
(ref. 24) in figures 2 to 4 (data shown as the solid curve) and in table 3.3.1. These functions 
should be improved in future work. 

3.4. Fragmentation Cross Sections 

The nuclear cross section for producing fragment Ay, Zj from a target At, Zt is taken from 
Rudstam (ref. 25) as 




l-5\ 

Z f -SA f +VA 2 f J 


D 


(3.4.1) 


where 


D = 1.79 

e vM ( 

j 0.3 \ 0.3 0.3' 

pA t ) A t pA t _ 

(3.4.2) 

Fi = U.8 2/3 e~ G+HAt 

(3.4.3) 

\e K 

f 2 = \ 

l io 

-LE 

(E < 240 MeV) 1 
(E > 240 MeV) J 

(3.4.4) 

( 20£-- 77 
P ~ { 0.056 

(E < 2100 MeV) 1 
(E > 2100 MeV) J 

(3.4.5) 

r = 11.8 AJ 0A6 


(3.4.6) 

S = 0.486 



(3.4.7) 

V = 3.8 x 10~ 4 


(3.4.8) 

G = 0.25 



(3.4.9) 
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H = 0.0074 


(3.4.10) 

(3.4.11) 

(3.4.12) 


K = 1.73 
L = 0.0071 

Correction factors for 12 C and 16 0 fragments are given in table 3.4.1. The total of all 
fragmentations is renormalized to the total absorption cross section at each energy. The major 
isotope fragments are reasonably well produced for light to medium-heavy fragments. Many 
other corrections to the basic Rudstam formalism have been made by Silberberg and Tsao 
(ref. 26). 

3.5. Differential Nuclear Cross Sections 

3.5.1. Nucleon-nucleon spectrum. The nucleon-nucleon differential cross sections are repre- 
sented (ref. 27) by 


f(E, E') = B + e ~ BE ] / (l - e~ BE ') (3.5.1) 

where 

B = 2mc 2 6/10 6 (3.5.2) 

In the above, me 2 is the nucleon rest energy (938 MeV) and b is the usual slope parameter, 
given by (units of GeV -2 ) 


( 3 + 14e- £ '/200 (f or ) ^ 

b={ , / (3.5.3) 

l 3.5 -I- 30e E / 200 (for np) J 

where E' is the initial nucleon energy in the rest frame of the target. The differential spectrum 
is defined (nonrelativistically) over the energy interval 0 < E < E' . Note that the expression 
(3.5.1) reduces to the usual result for low energy scattering 

f{E,E')*l/E' (3.5.4) 


The forward to backward scattering ratio is required for neutron scattering and is given by 
(ref. 28) 

F B {E') = 0.12 - .015# + - - J^-2.1) < 3 - 5 ' 5 ) 

where E' in equation (3.5.5) is the laboratory energy in GeV before collision. 

The differential cross sections are normalized such that 

j^=<r(E')f(E,E') (3.5.6) 

where a(E f ) is the appropriate nucleon-nucleon total cross section. Obviously we have neglected 
the inelastic processes, which must yet be included. The center of mass angular distributions 
{Ocm) are related to the energy change in the laboratory system by 

dn 4tt dE 

and are compared with the compilation of experimental data (ref. 29) in figures 5 and 6. These 
comparisons indicate that the present functions are reasonable. 
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3.5.2. Nucleon-nucleus spectrum. The nucleon-nucleus differential cross section in the 
impulse (Chew) approximation (note, this is just the Bom term of the optical model, ref. 30) 
is given by 

— = ce- 2b ^\F A (q 2 )\ 2 « ce -2b? 2 e -2a 2 g 2 /3 (3.5.8) 

where b is the slope parameter of equation (3.5.2) averaged among nuclear constituents, q is 
the magnitude of momentum transfer, and a is the nuclear rms radius. The nuclear rms radius 
in terms of the rms charge radius in fm is given as 


a = \J a* — 0.64 

where the rms charge radius in fm is 


r 0.84 

(At = 1) 

2.17 

II 

to 

1.78 

CO 

II 

< 1.63 

c-t- 

II 

4^ 

2.4 

(6 < A t < 14) 

, 0.82Al /3 + 0.58 

(At > 16) 


(3.5.9) 


(3.5.10) 


The nuclear form factor is the Fourier transform of the nuclear matter distribution. Note that 
the above assumes the nuclear matter distribution is a Gaussian function, which is reasonable 
for the light mass nuclei but is less valid for At » 20. 

The energy transferred to the nucleus Et is restricted by kinematics to 


where 


0 < E t < (1 - ot)Ef 
(At ~ l) 2 


a = 


(At + 1) 2 

The energy transfer spectrum is given as 

4Atmc 2 (B + a 2/3) e -4Atmc 2 (B+a 2 /3).Et 


f(E t ,E') = 


1 _ e -4Atmc 2 (l—a)(B+a?/3)E' 


(3.5.11) 

(3.5.12) 


(3.5.13) 


Similarly, the nucleon energy after scattering E is restricted to 

aE' <E< E' 


The nucleon spectrum is given by 


f(E,E') 


4A t mc 2 (B + a 2/ 3 ) e -4A t mc 2 (B+a 2 /3)(£'-.E) 
~ 1 _ e -4A t mc 2 (l-a)(B+a 2 /3)£'' 


(3.5.14) 


(3.5.15) 


One should note that both equations (3.5.13) and (3.5.15) reduce to the usual isotropic 
scattering result at low incident energy. 

The results of equation (3.5.15) are compared with experiment (refs. 31 and 32) in figures 7 
to 10. The comparison is rather good at the small angles, considering the simplicity of the 
present results, but there is definite room for improvement. Much of the present discrepancy 
is due to errors in o g (E) to which the present spectra are normalized in equation (3.5.7). 
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3.5.3. Nucleon nonelastic spectrum. The nonelastic differential cross sections use the results 
of Bertini’s MECC7 program. The nucleon multiplicities are given in tables 3.5.1 and 3.5.2. 
We have required the multiplicities to be monotonic in energy so the values in parentheses 
are corrected values obtained by scaling from lower and higher energies and are used in the 
calculations. The results below 400 MeV were taken from Alsmiller et al. (ref. 33), and the 
results for carbon, calcium, bromine cesium, and holmium above 400 MeV are obtained by 
interpolation. The nonelastic spectra are represented as 


f(E,E') 


y, Nj e~ g / a » 
iZ 1 1 ~ e~ E 'l a i 


+ 


E' 


Nq 

1 + e -20(l -E/E') 


(3.5.16) 


The first term of the summation represents the evaporation peak so that N\ is taken from 
table 3.5.1 and the spectral parameter ai (GeV) is taken from Ranft (ref. 34) as 


<*l p 


“in 


(0.019 + 0.0017£')(1 - 0.001 A t ) 
0.027(1 - O.OOlAf) 

(0.017 + 0.0017£')(1 - O.OOlAf) 
0.023(1 - O.OOlAt) 


(E' < 5 GeV) 1 
{E 1 > 5 GeV) J 

(E' < 5 GeV) 1 
(E' > 5 GeV) J 


(3.5.17) 


(3.5.18) 


The second term is taken from Ranft (ref. 34) to represent the low energy cascade particles as 


n 2p 


0.0035v/AI 

0.007v^ [0.5 + (1 + logo E ') 2 
0.0245 y/T t 


(E 1 < 0.1 GeV) 
(0.1 < E' < 5 GeV) 
( E ' > 5) 


(3.5.19) 


n 2 n 


r 0.0042\Mt 

< 0.007\/Xt [0.6 + 1.3(1 + log 10 E') 2 
k 0.032v^ 


with the corresponding spectral parameters 

_ / (0.11 + 0.01E') (1 - O.OOlAt) 
a2 P~ { 0.16(1 -O.OOlAt) 

_ / (0.1 + 0.01£') (1 - O.OOlAt) 
a2n ~ { 0.15(1 - O.OOlAt) 


(E' < 0.1 GeV) 1 
(0.1 < E' < 5 GeV) 1 


(3.5.20) 


(E' > 5 GeV) J 


(E' <5 GeV) 1 
[E' > 5 GeV) 

(E' < 5 GeV) 
(E 1 > 5 GeV) . 


(3.5.21) 


The third term in the summation is the balance of cascade particles after inclusion of the 
quasi-elastic contribution. 

The quasi-elastic contribution is estimated and includes the nuclear attenuation following 
the quasi-elastic event. The proton quasi-elastic cross section is 

VQ^pp = Zt&pp + {At — Zt)cr n p 

a Q,pn = {At — Z t )(j n p (3.5.22) 

and similarly for neutrons 

a Q,nn ~ (At Zt)(J nn + Zt&np 

^QyTip = Zt^np (3.5.23) 
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The corresponding multiplicities are taken as 

N Q,jk = e~' m ^ X *<TQ t jk / J2 a Qjl (3.5.24) 

t 

where the exponential factor accounts for the attenuation of the quasi-elastic particles before 
they escape the nucleus. The balance of the cascade particles are contained in N 3 as 


n 3 = n c -n 2 - n q 


(3.5.25) 


with an assumed spectral coefficient given by 

a 3 = a 2 /0.7 (3.5.26) 

Results of the present formalism are shown in figures 11 to 24 in comparison with the 
calculations of Bertini. Some further improvements in this parameterization need to be made. 


3.5.4 . Light fragment spectrum. The light fragment yields per event are given in table 3.5.3 as 
obtained from Bertini’s MECC7 calculations. These results are extrapolated and interpolated 
in energy and mass number. The corresponding mean energies are given in table 3.5.4. The 
mean energies are used in Ranft’s formula for nucleons and similarly for the light ions. 


3.5.5. Heavy fragment spectrum. The differential spectrum of the target fragments is obtained 
from the momentum distributions of projectile fragments produced in collision with hydrogen 
targets. The fragment densities in phase space were measured by the Heckman group (ref. 35) 
at Lawrence Berkeley Laboratory (LBL), from which the differential energy spectrum is found 
to be 

do o 0 exp{-E/2E 0 )\/E /0 c 0(rt 

which is normalized such that 

roc 

°0 = ~JB dE (3.5.28) 

JO dE 

is the fragmentation cross section. The value of Eq depends on the initial target mass and the 
fragment mass as 


£o 

b 


&A 


A p 


Ap/2 A f 



(3.5.29) 

min j 

(H2 A\ /2 

,260} 

(3.5.30) 

r 0.45 

(At = A f ) | 

(3.5.31) 


-A f 

(Otherwise) j 


' 46,4 

-j 

1/2 


0.86 

[20{A t - 

i). 

(MeV/c) 

(3.5.32) 


Values of A p obtained from the above are compared with experiments and results of others 
(ref. 35) in table 3.5.5. 

In relating the ultimate effects of the interaction of nucleons in materials, the energy 
given over to various ion fragment types per collision event is of utmost importance. The 
fragmentation and energy transfer cross sections of Bertini and the present calculation are 
compared with the experiments of Greiner et al. in table 3.5.6. The energy transfer of the 
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Bertini data set is nearly half of that observed in experiments. The present result with the 
Rudstam cross sections is low by less than 20 percent. An ad hoc correction factor is applied 
to the Rudstam cross sections as shown in table 3.5.6 to better match the experimental data. 
Note, the corrections make the results slightly conservative. 

4. Results 

In an effort to begin validation of the present code, we compare the results with prior 
calculations using Monte Carlo methods and essentially the same data base. Fully three- 
dimensional Monte Carlo calculations have been made with the Bertini code as the nuclear 
cross section set and low energy neutron data. (See refs. 36 and 37 for a detailed discussion.) 
Energy absorption in a tissue slab for normally incident neutrons of energies 0.5, 2, and 10 MeV 
is shown in table 4.1. Also shown are the results of the present code. The results appear 
remarkably good considering the crudeness of the straight ahead approximation for low energy 
neutrons and the limitations of the present data base. Results for higher energy neutrons are 
shown in figures 25 to 29. In each case, reasonable agreement with the results of Zerby and 
Kinney is obtained. Similar results are found for energetic protons as shown in figures 30 
to 33. In the present calculation, the first generation proton spectrum is discontinuous for 
monoenergetic beams and is best handled by taking many energy points in the spectrum. 
However, the calculation time then becomes excessive. The present results were calculated 
using only 30 energy points; this is adequate for space radiation, as shown in reference 38, 
but marginal for the present monoenergetic results. The use of numerical benchmark problems 
will allow us to better understand the numerical procedures. Such a benchmark has already 
provided some insight (ref. 38). 

5. Concluding Remarks 

The emphasis of the present report is on high energy baryon transport and a relatively 
complete model is presented. But it is especially important that any such code adequately 
represent the low energy neutrons in a reasonable way. It is seen from these results that this 
has been accomplished in the present code employing the straight ahead approximation. The 
calculated doses of 100 to 400 MeV neutrons and protons on tissue are in reasonable agreement 
with the more complete Monte Carlo code. The primary advantage of the present code is 
computer efficiency while maintaining adequate accuracy. Future work will concentrate on 
improving the data base to further improve the comparisons. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
May 6, 1988 
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Table 3.3.1. Simplified Total Neutron-Nucleus Cross Sections and Fitted Data at 20 MeV 



Cross section, mb, for — 

6 Li 

9 Be 

n B 

160 

27 A1 

^Cu 

131 j 

209p b 

^exp 

Equation (3.3.5) 

1250 

1157 

1300 

1409 

1490 

1547 

1600 

1833 

1850 

2287 

2500 

3279 

5000 

4559 

5500 

5692 


Table 3.4.1. Ad Hoc Correction Factors for Rudstam’s Formula 


AA 

Correction factor for — 

12C 

16q 

1 

1.3 

1.5 

2 

0.5 

1.0 

3 

0.3 

1.0 

4 

0.1 

1.0 

5 

1.0 

1.5 

6 

0.35 

0.5 

7 

! 

0.5 

8 


0.5 

9 


0.5 

10 


1.0 
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Table 3.5.1. Number of Evaporation Nucleons Produced in Nuclear Collisions 
[Values in parentheses are modified and used in the code] 


Number of nucleons produced at- 




n — * p 
n — ► n 
A t = 16: 

P^P 
p — ► n 
n — > p 
n — ► n 
A t = 27: 

P^P 
p — ► n 
n — ► p 
n — > n 
A t = 40: 

p-+p 

p — ► n 
n — ► p 
n — * n 
= 65: 

P^P 
p — ► n 
n —+ p 
n — ► n 
A t = 80: 

P^P 
p n 
n — ► p 
n — ► n 
= 100 : 

P^P 
p — ► n 
n — > p 
n — * n 
At = 132: 

P^P 
p — * n 
n — ► p 
n — ► n 
At = 164: 

p->p 

p — ► n 
n — ► p 
n — ► n 
A t = 207: 

P-+P 
p — ► n 
n —> p 
n—>n 


25 MeV 

200 MeV 

400 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

0.51 

0.54 

0.50 

0.72 

0.75 

0.84 

0.026 

0.32 

0.35 

0.79 

0.79 

0.79 

0.052 

0.30 

0.35 

0.73 

0.73 

0.80 

0.43 

0.57 

0.52 

0.77 (0.71) 

0.71 (0.71) 

0.73 

0.62 

0.73 

0.71 

0.84 

0.89 

0.98 (0.93) 

0.07 

0.36 

0.441 

0.11 (0.87) 

0.93 (0.87) 

0.82 (0.87) 

0.12 

0.47 

0.53 

0.86 

0.86 

0.89 

0.55 

0.60 

0.59 

0.79 

0.79 

0.81 

0.54 

0.99 

1.03 

1.36 

1.49 

1.86 

0.37 

0.61 

0.62 

1.29 

2.03 (1.92) 

1.52 (1.92) 

0.14 

0.78 

0.82 

1.29 

1.60 

1.74 

0.75 

0.76 

0.71 

1.34 

... 1 51 

1.60 

0.50 

1.03 

1.06 

1.74 

2.32 

2.93 

0.53 

1.12 

1.24 

2.63 

3.36 

3.64 

0.12 

0.74 

0.84 

1.60 

2.29 

2.67 

0.89 

1.39 

1.44 

2.76 

3.25 

3.54 

0.18 

0.75 

0.91 

2.11 

3.15 

4.00 

1.04 

2.33 

2.65 

3.97 

4.79 

5.37 

0.03 

0.49 

0.66 

1.90 

2.98 

3.61 

1.46 

2.77 

2.90 

4.17 

4.99 

5.49 

0.10 

0.60 

1.07 

2.2 

3.18 

4.89 

1.29 

2.20 

3.18 

3.72 

5.07 

6.77 

0.02 

0.53 

0.79 

1.87 

2.91 

0.53 

1.58 

3.19 

3.43 

4.07 

5.35 

6.91 

0.03 

0.46 

1.28 

2.96 

4.56 

5.78 

1.53 

1.97 

3.72 

5.46 

7.04 

8.17 

0.004 

0.59 

0.96 

2.71 

4.27 

5.44 

1.67 

3.60 

3.97 

5.63 

7.31 

8.33 

0.01 

0.61 

1.03 

2.68 

4.51 

6.32 

1.91 

4.11 

5.25 

8.76 

11.34 

12.31 

0.001 

0.47 

0.81 

2.51 

4.47 

5.98 

1.96 

4.73 

5.59 

8.93 

10.6 

12.42 

0.003 

0.42 

0.76 

2.38 

4.68 

6.86 

2.17 

5.79 

7.07 

12.09 

15.7 

16.45 

0.003 

0.28 

0.58 

2.30 

4.68 

6.52 

2.26 

5.96 

7.07 

12.3 

14.6 

16.51 

0.001 

0.21 

0.44 

2.23 

5.19 

7.39 

2.29 

7.22 

9.24 

15.3 

17.81 

20.6 

0.00 

0.10 

0.30 

2.10 

4.88 

7.05 

2.29 

7.38 

9.53 

15.6 

18.2 

20.6 
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Table 3.5.2. Number of Cascade Nucleons Produced in Nuclear Collisions 


Number of nucleons produced at- 


A t = 12: 

P^P . 
p — ► n . 
n — * p . 
n n . 
A* = 16: 

. 

p —> n . 
n — ► p . 
n — > n . 
A t = 27: 

P^P . 
p — ► n . 
n — ► p . 
n — ► n . 
>17 = 40: 

p^p . 

p — ► n . 
n — ► p . 
n — ► n . 
^7 =65; 

p^p . 
p — ► n . 
n — ► p . 
n — ► n . 
i4 t = 80: 

p^p . 

p — ► n . 
n p . 
n — > n . 
= 100: 

P P . 
p — ► n . 
n p . 
n — ► n . 
^ =132: 

p^p . 

p — ► n . 
n — ► p . 
n — ► n . 
Ai = 164: 

p^p . 

p — ► n . 
n — ► p . 
n — ► n . 

: = 208: 
p-^p . 

p — ► n . 
n —> p . 
n — > n . 


25 MeV 

200 MeV 

400 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

0.58 

1.43 

1.63 

1.95 

2.15 

2.48 

0.41 

0.86 

0.93 

1.42 

1.66 

2.08 

0.42 

0.90 

0.92 

1.43 

1.65 

1.91 

0.56 

1.42 

1.69 

1.95 

2.27 

2.57 

0.56 

1.41 

1.72 

2.05 

2.39 

2.60 

0.38 

0.90 

0.98 

1.47 

1.86 

2.19 

0.38 

0.91 

0.96 

1.49 

1.85 

2.01 

0.54 

1.43 

1.70 

2.05 

2.52 

2.70 

0.46 

1.38 

1.67 

2.29 

2.86 

3.19 

0.34 

0.97 

1.16 

1.86 

2.54 

3.25 

0.32 

0.93 

1.01 

1.69 

2.28 

2.71 

0.49 

1.48 

1.81 

2.42 

3.22 

3.71 

0.40 

1.33 

1.69 

2.32 

3.01 

3.53 

0.30 

1.04 

1.24 

2.46 

3.52 

4.48 

0.28 

0.89 

1.08 

1.79 

2.51 

3.06 

0.45 

1.49 

1.88 

2.99 

4.13 

4.83 

0.30 

1.21 

1.69 

2.35 

3.16 

3.87 

0.28 

1.09 

1.46 

3.06 

4.49 

5.72 

0.21 

0.86 

1.08 

1.88 

2.75 

3.41 

0.40 

1.53 

2.00 

3.55 

5.03 

5.95 

0.27 

1.18 

1.57 

2.32 

3.18 

3.95 

0.25 

1.08 

1.45 

3.27 

4.92 

6.35 

0.19 

0.81 

1.04 

1.86 

2.78 

3.54 

0.36 

1.51 

1.98 

3.78 

5.40 

6.64 

0.25 

1.15 

1.55 

2.29 

3.20 

4.04 

0.22 

1.06 

1.52 

3.47 

5.35 

6.98 

0.17 

0.78 

1.08 

1.84 

2.44 

3.67 

0.31 

1.47 

2.03 

3.96 

5.76 

7.33 

0.20 

1.00 

1.46 

2.21 

3.17 

3.87 

0.20 

1.11 

1.57 

3.31 

5.20 

7.91 

0.13 

0.70 

1.00 

1.79 

2.69 

3.52 

0.28 

1.45 

2.10 

3.86 

6.86 

8.29 

0.16 

0.90 

1.36 

2.13 

3.15 

3.69 

0.18 

1.11 

1.60 

3.16 

5.06 


0.11 

0.63 

0.88 

1.72 

2.55 


0.26 

1.42 

2.11 

3.56 

7.94 


0.14 

0.82 

1.27 

2.05 

7.74 

3.51 

0.16 

1.03 

1.71 

2.97 

7.23 

9.77 

0.09 

0.58 

0.87 

1.67 

2.41 

3.24 

0.23 

1.36 

2.10 

3.36 

7.63 

10.21 


Table 3.5.3. Evaporated Ion Yields From Nucleon-Nucleus Collisions 
[Values in parentheses are for proton reactions] 




Ion yields at — 



500 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

A t = 16 : 





d 

0.111 (0.094) 

0.199 (0.237) 

0.257 (0.265) 

0.304 (0.311) 

t 

0.022 (0.029) 

0.024 (0.025) 

0.033 (0.025) 

0.029 (0.029) 

he 

0.018 (0.034) 

0.035 (0.043) 

0.037 (0.052) 

0.037 (0.048) 

a 

0.664 (0.400) 

0.720 (0.696) 

0.664 (0.624) 

0.640 (0.667) 

11 





d 

0.126 (0.130) 

0.245 (0.269) 

0.380 (0.396) 

0.442 (0.433) 

t 

0.028 (0.023) 

0.048 (0.052) 

0.063 (0.065) 

0.072 (0.069) 

he 

0.042 (0.035) 

0.067 (0.074) 

0.073 (0.091) 

0.083 (0.092) 

a ..... . 

0.370 (0.400) 

0.550 (0.566) 

0.597 (0.582) 

0.577 (0.577) 

A t = 65 : 





d 

0.150 (0.171) 

0.379 (0.390) 

0.748 (0.766) 

0.935 (0.987) 

t 

0.031 (0.035) 

0.075 (0.068) 

0.145 (0.145) 

0.177 (0.191) 

he 

0.013 (0.014) 

0.039 (0.056) 

0.112 (0.124) 

0.166 (0.177) 

a 

0.124 (0.137) 

0.231 (0.231) 

0.373 (0.377) 

0.431 (0.441) 

A t = 100 : 





d 

0.174 (0.183) 

0.456 (0.475) 

1.01 (1.02) 

1.44 (1.48) 

t 

0.028 (0.029) 

0.080 (0.081) 

0.207 (0.192) 

0.269 (0.273) 

he 

0.012 (0.017) 

0.055 (0.060) 

0.162 (0.185) 

0.249 (0.262) 

a 

0.158 (0.156) 

0.320 (0.339) 

0.490 (0.467) 

0.549 (0.540) 

A t = 207 : 
d 

0.131 (0.i52) 

0.536 (0.565) 

1.51 (1.57) 

2.54 (2.54) 

t 

0.038 (0.037) 

0.152 (0.163) 

0.415 (0.424) 

0.641 (0.644) 

he 

0.001 (0.002) 

0.017 (0.017) 

0.112 (0.106) 

0.211 (0.239) 

a 

0.053 (0.063) 

0.195 (0.210) 

0.527 (0.514) 

0.751 (0.746) 
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Table 3.5.4. Mean Energies (MeV) of Light Nuclear Fragments Produced in Nucleon-Nucleus Collisions 

[Values in parentheses are for proton reactions] 



Mean energies at — 


500 MeV 

1000 MeV 

2000 MeV 

3000 MeV 

A t = 16: 





n 

5.55 (6.19) 

7.91 (7.89) 

9.55 (9.81) 

11.1 (9.80) 

P 

6.10 (6.40) 

8.33 (8.69) 

9.71 (10.2) 

10.3 (11.2) 

d 

8.53 (7.64) 

12.2 (10.7) 

14.9 (14.8) 

16.3 (13.0) 

t 

6.40 (7.83) 

10.6 (10.4) 

12.5 (9.74) 

13.7 (10.1) 

he 

12.1 (8.76) 

11.8 (11.2) 

11.1 (13.1) 

12.9 (10.3) 

a 

9.36 (6.24) 

12.6 (12.3) 

13.1 (14.6) 

13.6 (13.8) 

A t = 27: 





n 



9.91 (10.5) 


P 



11.1 (11.9) 

13.5 (13.7) 

d 

9.57 (9.42) 


14.3 (14.8) 

17.2 (17.4) 

t 

9.16 (9.54) 


13.0 (13.9) 

16.6 (13.7) 

he 



13.4 (14.1) 

14.4 (14.5) 

a 

12.7 (13.4) 


13.8 (13.8) 

14.5 (14.6) 

A t = 65: 





n 

4.24 (4.32) 

5.67 (5.70) 

7.92 (7.91) 

9.67 (9.58) 

p 

8.25 (8.30) 

9.66 (9.76) 

12.1 (12.3) 

14.4 (14.2) 

d 

9.88 (10.1) 

13.5 (11.8) 

13.8 (14.2) 

15.6 (15.9) 

t 

10.0 (10.0) 

11.7 (11.6) 

13.7 (13.8) 

15.1 (15.9) 

he 

14.6 (14.1) 

16.4 (16.2) 

17.5 (19.3) 

19.5 (19.2) 

a 

12.7 (13.4) 

13.2 (13.6) 

13.8 (13.8) 

14.5 (14.6) 

o 

o 

i-H 

II 





n 

3.90 (3.90) 

5.13 (5.16) 

7.11 (7.04) 

8.61 (8.74) 

p 

9.63 (9.62) 

11.0 (11.0) 

12.9 (13.2) 

14.6 (14.7) 

d 

11.0 (11.1) 

12.5 (12.6) 

14.4 (15.0) 

16.1 (16.0) 

t 

11.3 (11.7) 

12.6 (13.0) 

14.7 (14.3) 

15.5 (16.5) 

he 

17.8 (18.7) 

18.6 (18.8) 

20.9 (20.6) 

21.8 (22.2) 

a 

16.5 (16.5) 

16.8 (16.9) 

17.5 (17.5) 

17.6 (17.6) 






n 

3.28 (3.27) 

4.37 (4.33) 

5.83 (5.78) 

6.90 (6.95) 

p 

12.5 (12.5) 

12.2 (13.4) 

14.9 (14.9) 

16.2 (16.3) 

d 


14.4 (14.2) 

16.0 (16.8) 

17.4 (17.8) 

t 


5.0 (15.3) 

16.6 (16.8) 

17.4 (17.8) 

he 


26.2 (26.5) 

28.0 (27.8) 

29.1 (28.5) 

a 

25.3 (25.7) 


26.4 (26.3) 

25.9 (26.4) 
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Table 3.5.5. A p for 16 0 Fragments Produced by 2.1 GeV Protons 


Fragment 

Value of A p , MeV/c, from — 

Experiments (ref. 35) 

Present work 

Greiner 

rsQ 

94 ±3 

80.0 

83.8 

14 o 

99 ±6 

109.5 

113.1 

13q 

143 ± 14 

129.2 

133.5 

16 n 

54 ±11 

55.0 

? 

15 n 

95 ±3 

80.0 

82.8 

14 n 

112 ±3 

109.5 

113.0 

13 n 

134 ±2 

129.2 

133.5 

12 n 

153 ±11 

143.4 

148.1 

i5 C 

125 ± 19 

80.0 

82.8 

14 c 

125 ±3 

109.5 

113.1 

13 c 

130 ±3 

129.2 

133.5 

12 C 

120 ±4 

143.36 

148.09 

n c 

162 ±5 

153.45 

158.5 

10 c 

190 ±9 

160.3 

165.6 

13 b 

166 ± 10 

129.2 

133.5 

12b 

163 ±8 

143.4 

148.1 

n B 

160 ±2 

153.5 

158.5 

i°B 

175 ±7 

160.3 

165.6 

8 B 

175 ± 22 

165.5 

171.01 

“Be 

197 ± 20 

153.5 

158.5 

1° Be 

159 ±6 

160.0 

165.0 

7 Be 

166 ±2 

164.24 

169.66 

9 Li 

188 ±15 

164.24 

169.66 

8 Li 

170 ± 13 

165.4 

171.0 

7 Li 

163 ±4 

164.2 

169.66 

6 He 

167 ± 20 

1 

160.0 

165.0 
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Table 3.5.6. Fragmentation Cross Sections and Fragment Energy Transfer Cross Sections 
From Rudstam, Bertini, and Experiments (Grainer) 

[Values in parentheses include the ad hoc factor] 


a f 

Fragmentation cross sections, 
cr, mb, from — 

Energ 

J 

ry transfer cross sections, 
5(7, MeV-mb, from — 

Ad hoc 
factor 

— 

Bertini 

Greiner 

Rudstam 

Bertini 

Greiner 

Rudstam 

16 

4.69 



5.04 

0.0006 



15 

103.4 

61.5 

24.8 

60.6 

56.9 

17.0 


14 

40.0 

35.4 

23.6 

48.8 

51.7 

32.3 


13 

18.5 

22.8 

22.4 

37.6 

48.3 

46.0 


12 

32.2 

34.1 

22.0 

85.8 

68.2 

60.3 


11 

8.2 

26.4 

17.1 

37.9 

99.1 

58.2 


10 


12.7 

20.6 

52.8 

62.0 

84.6 


9 

1.2 

5.2 

12.8 

6.5 

25.7 

61.4 

0.5 

8 

.56 

1.23 

1.87 

2.5 

7.1 

10.4 


7 

.16 

27.9 

9.1 

6.11 

153.4 

56.1 

3 

6 

5.46 

17.5 

18.0 

31.4 

92.5 

123.3 


Total 

225.4 

244.8 

172.3 (265.8) 

375.1 

664.9 

549.6 (741) 
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Table 4.1. Energy Deposition of 0.5-10 MeV Neutrons 
[Values in parentheses are from present calculations] 




Energy deposition from Monte Carlo and present calculations, MeV 

Incident 






energy, MeV 

Depth, cm 


Proton 


Heavy ion 

0.5 

z <*> 

0.1107 

.0986 

0.0767 

0.0073 

.0066 

0.0052 


«•> 

.0418 

.0331 

(.0296) 

.0028 

.0018 

(.0021) 


^ <»» 

.0074 

.0059 

.0085) 

.0003 

.0002 

(.0005) 


14- 15 , , 

15- 16 ^ 

.0006 

.0007 

(.0016) 


(.0001) 


r <“> 

.0002 

.0001 

(.0005) 



2 

z <*> 

0.2138 

.1984 

(0.1773) 

0.0147 

.0133 

(0.0124) 


s «•> 

.1539 

.1349 

(.1122) 

.0110 

.0105 

(.0074) 


M0 (.0) 
10-11 v } 

.0770 

.0741 

(.0602) 

.0054 

.0061 

(.0045) 


14 ~ 15 (15) 

15-16 } 

.0301 

.0240 

(.0301) 

.0026 

.0019 

(.0023) 


1M0 (20) 
20-21 V ' 

.0091 

.0114 

(.0142) 

.0008 

.0009 

(.0011) 

10 

Z <» 

0.3520 

.3284 

(0.3056) 

0.0339 

.0342 

(0.0376) 


n <•> 

.3220 

.2977 

(.2249) 

.0345 

.0266 

(.0267) 


<»» 

.2674 

.2661 

(.1537) 

.0257 

.0225 

(.0176) 


:ti: <»> 

.2161 

.2211 

(.1049) 

.0203 

.0198 

(.0108) 


™ <*» 

.1635 

.1291 

(.0714) 

.0150 

.0149 

(.0078) 
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BichseKref. 15) 

Palmerand Akhavan-Rezayat (ref. 20) 
Matteson et al.(ref. 19) 

Northcliffe and Schilling(ref. 21) 
Calculated 



o/Aei/M ‘J9Mod 6u|ddois 


kinetic energy. 





E, MeV 

Figure 2. Total cross section of neutrons with lithium-6 as represented in the present code and according to 
experiment (ref. 24). 



0.1 1 10 100 


E, MeV 

Figure 3. Total cross section for neutrons with carbon-12 according to present code and experiments (ref. 24). 
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E, MeV 


Figure 4. Total cross section for neutrons with aluminum-27 according to present code and experiments 
(ref. 24). 



Figure 5. Neutron-proton scattering cross sections according to present code and experiment. 




0 ° 30 ° 60 ° 90 ° 

0 

cm 

Figure 6. Proton-proton scattering cross sections according to present code and experiment. 



Cos 0 


cm 


Figure 7. Neutron-carbon scattering cross sections according to present code and experiment (ref. 32). 
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238 

n- U 
BNL 400 


Present 


Figure 8. Neutron-uranium scattering cross sections according to present code and experiment (ref. 32). 


d a 
dO ’ 

millibarns/sr 



7 MeV Neutrons 

Fernbach 

Present 


o 


Figure 9. Neutron scattering cross sections for several elements according to present code and experiments. 
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Figure 10. Neutron scattering cross sections for several elements according to present code and experiments. 



Figure 11. Nucleon cascade spectra of Bertini and the present work. Protons produced by 100 MeV protons 
on oxygen. 
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Figure 12. Nucleon cascade spectra of Bertini and the present work. Protons produced by 500 MeV protons 
on oxygen. 



Figure 13. Nucleon cascade spectra of Bertini and the present work. Protons produced by 1000 MeV protons 
on oxygen. 



Figure 14. Nucleon cascade spectra of Bertini and the present work. Protons produced by 3000 MeV protons 
on oxygen. 



Figure 15. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 1000 MeV protons 
on oxygen. 
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Secondary neutron energy E, MeV 


Figure 16. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 3000 MeV protons 
on oxygen. 



Figure 17. Nucleon cascade spectra of Bertini and the present work. Protons produced by 1000 MeV protons 
on aluminum. 





Figure 18. Nucleon cascade spectra of Bertini and the present work. Protons produced by 3000 MeV protons 
on aluminum. 



Figure 19. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 1000 MeV protons 
on aluminum. 
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Present 


Secondary neutron energy E, MeV 

Figure 20. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 3000 MeV protons 
on aluminum. 


f(E,E 0 ). 



Secondary proton energy E, MeV 


Figure 21. Nucleon cascade spectra of Bertini and the present work. Protons produced by 1000 MeV protons 
on lead. 






Secondary proton energy E, MeV 


Figure 22. Nucleon cascade spectra of Bertini and the present work. Protons produced by 3000 MeV protons 
on lead. 



Figure 23. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 1000 MeV protons 
on lead. 




Figure 24. Nucleon cascade spectra of Bertini and the present work. Neutrons produced by 3000 MeV protons 
on lead. 


Dose, 

rad/neutron/cm 



Figure 25. Dose in tissue due to normal incident 60 MeV neutrons according to reference 37 and present code 





Depth, cm 


Figure 26. Dose in tissue due to normal incident 100 MeV neutrons according to reference 37 and present code. 


Dose, 

rad/neutron/cm 



Depth, cm 


Figure 27. Dose in tissue due to normal incident 200 MeV neutrons according to reference 37 and present code. 
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Depth, cm 


Figure 28. Dose in tissue due to normal incident 300 MeV neutrons according to reference 37 and present code. 


Dose, 

rad/neutron/cm 



Figure 29. Dose in tissue due to normal incident 400 MeV neutrons according to reference 37 and present code. 
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10 


6 



Depth, cm 


Figure 30. Dose in tissue due to normal incident 100 MeV protons according to reference 37 and present code. 
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Dose, 

rad/proton/cm 



Figure 31. Dose in tissue due to normal incident 200 Mev protons according to reference 37 and present code. 
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Figure 32. Dose in tissue due to normal incident 300 MeV protons according to reference 37 and present code. 



Figure 33. Dose in tissue due to normal incident 400 MeV protons according to reference 37 and present code. 
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